NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 

EFFECTS  OF  WATER  VAPOR  AND  ANISOTROPIC 

SCATTERING  ON  AEROSOL 

OPTICAL  DEPTH  RETRIEVAL 

by 

Mary  B.  Clifford 

December  1991 

Thesis 

Advisor                     Philip  A. 

Du] 

:kee 

Approved  for  public  release;  distribution  unlimited 


T257794 


nclassified 


cunty  classification  of  this  page 


REPORT  DOCUMENTATION   PACE 

a  Report  Security  Classification  Unclassified 

lb  Restrictive  Markings 

a  Security  Classification  Authority 

3  Distribution  Availability  of  Report 

Approved  for  public  release;  distribution  is  unlimited. 

b  Declassification  Downgrading  Schedule 

Performing  Organization  Report  Number(s) 

5  Monitoring  Organization  Report  Number's) 

a  Name  of  Performing  Organization 

<aval  Postgraduate  School 

6b  Office  Symbol 
(if  applicable)  35 

7a  Name  of  Monitoring  Organization 

Naval  Postgraduate  School 

c  Address  (dry,  state,  and  ZIP  code) 

Monterey.  CA  93943-5000 

7b  Address  (city,  state,  and  ZIP  code) 

Monterey.  CA  93943-5000 

a  Name  of  Funding  Sponsoring  Organization 

8b  Office  Symbol 
(if  applicable) 

9  Procurement  Instrument  Identification  Number 

c  Address  (city,  state,  and  ZIP  code) 

10  Source  of  Funding  Numbers 

Program  Element  No 

Project  No 

Task  No 

Work  Unit  Accession  No 

1  Title  (include  security  classification) 

EFFECTS  OF  WATER  VAPOR  AND  ANISOTROPIC  SCATTERING  ON  AEROSOL 

OPTICAL  DEPTH  RETRIEVAL 

2  Personal  Author(s)  Man'  B.  Clifford 

3a  Type  of  Report 
Master's  Thesis 

13b  Time  Covered 
From                        To 

14  Date  of  Report  (year,  month,  dav) 

December  1991 

15  Page  Count 
48 

6  Supplementary  Notation  The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official  policy  or  po- 
rtion of  the  Department  of  Defense  or  the  U.S.  Government. 

7  Cosati  Codes 

18  Subject  Terms  (continue  on  reverse  if  necessary  and  identify  by  block  number) 

word  processing,  Script,  GML,  text  processing. 

ield 

Group 

Subgroup 

9  Abstract  (continue  on  reverse  if  necessary  and  identify  by  block  number) 

NOAA-7  AVIIRR  data  from  April  1982  was  used  to  perform  global-scale  analysis  of  aerosol  particle  characteristics, 
•-lahony's  (1991)  water  vapor  correction  was  incorporated  into  the  AVE1RR  multichannel  satellite  data  processing  technique 
sed  by  Frost  (1988).  Channel  4  Channel  5  brightness  temperature  difference  was  used  as  an  estimate  of  water  vapor  content 
i  an  air  column.  Greatest  measured  water  vapor  content  was  at  the  equator,  decreasing  toward  the  poles.  Applying  the  cor- 
ection  reduced  the  aerosol  particle  size  index,  resulting  in  an  increase  in  aerosol  optical  depths.  In  addition,  the  high  particle 
ize  index  over  the  low  latitudes  in  the  southern  hemisphere  noted  by  Frost,  was  reduced  significantly  after  applying  the 
orrection.  A  comparison  was  made  between  Frost's  values  for  the  variable  two-term  Henyey-Greenstein  phase  function  and 
lose  derived  with  the  water  vapor  correction.  The  basic  shape  of  the  curve  is  similar  to  that  determined  by  Frost. 

0  Distribution  Availability  of  Abstract 

3  unclassified  unlimited          □  same  as  report          □  DTIC  users 

21  Abstract  Security  Classification 

Unclassified 

2a  Name  of  Responsible  Individual 
'hilip  A.  Durkee 

22b  Telephone  (include  Area  code) 
(408)  646-3465 

22c  Office  Symbol 
MR  De 

D  FORM  1473.84  MAR 


83  APR  edition  may  be  used  until  exhausted 
All  other  editions  are  obsolete 


security  classification  of  this  page 


Approved  for  public  release;  distribution  is  unlimited. 

Effects  of  Water  Vapor  and  Anisotropic  Scattering  on  Aerosol 
Optical  Depth  Retrieval 


bv 


Mary  B.  Clifford 

Lieutenant,  United  States  Navy 

B.S.,  State  Univeristy  of  New  York  at  Albany,  1979 

Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  METEOROLOGY  AND  PHYSICAL 

OCEANOGRAPHY 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 

December  1991 


,    A 


-i 3?* — —      r 

Robert  L.  Haney,  Chairma^f 
Department  of  Meteorology 


ABSTRACT 

NOAA-7  AVHRR  data  from  April  1982  was  used  to  perform  global-scale  analysis 
of  aerosol  particle  characteristics.  Mahony's  (1991)  water  vapor  correction  was  incor- 
porated into  the  AVHRR  multichannel  satellite  data  processing  technique  used  by  Frost 
(1988).  Channel  4/Channel  5  brightness  temperature  difference  was  used  as  an  estimate 
of  water  vapor  content  in  an  air  column.  Greatest  measured  water  vapor  content  was 
at  the  equator,  decreasing  toward  the  poles.  Applying  the  correction  reduced  the  aerosol 
particle  size  index,  resulting  in  an  increase  in  aerosol  optical  depths.  In  addition,  the  high 
particle  size  index  over  the  low  latitudes  in  the  southern  hemisphere  noted  by  Frost,  was 
reduced  significantly  after  applying  the  correction.  A  comparison  was  made  between 
Frost's  values  for  the  variable  two-term  Henyey-Greenstein  phase  function  and  those 
derived  with  the  water  vapor  correction.  The  basic  shape  of  the  curve  is  similar  to  that 
determined  bv  Frost. 
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I.     INTRODUCTION 

A.  PURPOSE 

The  availability  of  satellite-measured  radiances  provides  an  effective  and  relatively 
inexpensive  means  for  determining  aerosol  particle  distribution  on  a  global  scale.  Satel- 
lites provide  extensive  coverage  of  the  earth's  surface  with  sufficient  temporal  periodicity 
to  allow  for  the  study  of  large-scale  global  effects  while  at  the  same  time  providing 
small-scale  local  measurements  desirable  for  military  applications.  Data  from 
meteorological  satellites  require  careful  processing  in  order  to  distinguish  the  much  lower 
radiances  associated  with  aerosol  optical  depth  variations  (Durkee  1984)  from  the 
greater  radiances  associated  with  cloud  albedos,  satellites  were  designed  to  measure. 

Improved  aerosol  optical  depth  retrieval  algorithms  have  facilitated  ongoing  efforts 
to  establish  a  global  aerosol  particle  data  base.  Frost  (1988)  developed  an  algorithm  to 
derive  aerosol  particle  characteristics  from  the  Advanced  Very  High  Resolution 
Radiometer  (AVHRR)  on  NOAA  polar  orbiting  satellites.  Aerosol  particle  size  index 
can  be  derived  by  a  comparison  of  the  sensitivities  of  two  adjacent  spectral  windows  and 
so  detect  aerosol  scattering.  Frost  examined  the  performance  of  a  variable  phase  func- 
tion as  an  improvement  to  the  constant  phase  function  used  previously.  Mahony  (1991) 
examined  the  effects  of  water  vapor  on  aerosol  optical  depth  retrieval.  An  estimate  of 
water  vapor  content  in  a  column  of  air  can  be  derived  by  comparing  the  brightness 
temperature  of  two  AVHRR  windows  in  the  thermal  infrared  region.  He  concludes  that 
this  effect  can  be  compensated  for,  eliminating  an  error  of  15%  in  the  phase  function. 

The  purpose  of  this  study  is  to  incorporate  Mahony's  improvement  into  the  algo- 
rithm used  by  Frost,  which  did  not  account  for  water  vapor  appreciably  affecting  re- 
trieved radiances.  The  April  19S2  data  examined  by  Frost  will  be  re-evaluated  to 
compensate  for  water  vapor.  The  variable  two-term  Henyey-Greenstcin  phase  function 
recommended  by  Frost  will  be  examined  in  light  of  the  new,  water  vapor  corrected  par- 
ticle size  phase  function  and  compared  to  the  values  he  generated. 

B.  MOTIVATION 

Aerosols  consist  of  a  variety  of  solid  or  liquid  material  dispersed  in  the  atmosphere; 
examples  include  volcanic  ash,  smoke,  and  pollen.  These  particles  vary  in  their  distrib- 
ution and  concentration  in  the  atmosphere,  which  can  impact  local  and  global  radiation 
budgets.    The  extreme  concentrations  of  smoke  particles  resulting  from  the  oil  fires  in 


Kuwait  created  a  layer  that  reduced  daylight  to  night  in  the  region.  Even  modest  con- 
centrations have  demonstrated  an  impact  in  models  of  the  earth's  radiation  budget 
(Weiss  et  al.  1974;  and  Coakley  1976).  A  Soviet  study  of  longwave  radiation  in  the 
tropical  Atlantic  ocean  (Zaitseva  1976)  linked  short-term  temporal  changes  in  temper- 
ature to  changes  in  aerosol  particles. 

Charlson  et  al.  (1987)  proposed  that  the  major  source  of  cloud  condensation  nuclei 
(CCN)  over  remote  unpolluted  oceans  appears  to  be  a  waste  product  of  phytoplankton, 
aqueous  dimethylsulphide  gas  (DMS).  A  complex  chain  of  chemical  and  biological 
processes  result  in  DMS  passing  through  the  sea  surface.  Once  in  the  atmosphere,  it  is 
oxidized  to  form  sulfate  aerosol  particles  which  act  as  CCN.  Varitions  in  phytoplankton 
populations  with  the  associated  changes  in  CCN  particles  could  alter  the  earth's  heat 
budget. 

Of  direct  importance  to  the  military  is  the  impact  of  aerosols  on  electro-optical  de- 
vices. The  .4  to  1.0  (im  waveband  is  considered  the  optimum  window  for  lasers  as  it 
avoids  ozone  absorption  at  shorter  wavelengths  and  water  vapor  absorption  at  longer 
wavelengths  (Blocmbergen  et  al.  1987).  Aerosol  scattering  is  a  significant  factor  in  this 
same  waveband. 


II.     THEORY 

A.     DIFFUSE  RADIANCE 

Knowing  that  optical  depths  are  small  over  the  ocean  and  that  single  scattering  is 
dominant  (Pfeil  19S6),  and  assuming  an  optically  thin  atmosphere.  Frost  (19S8)  used  the 
following  simplified  version  of  the  Radiation  Transfer  Equation  (RTE): 


L{5,ix,$)~(^^X{®)5,  2.1 


where, 

L  =  diffuse  radiance, 

4>  =  azimuth  angle, 

<y0  =  single  scatter  albedo, 

F0  =  solar  irradiance, 

H  =    cos  6,  (6  =  satellite  zenith  angle), 

p  =  scattering  phase  function, 

0  =  single  scatter  angle, 

S  =  optical  depth. 

Griggs  (1975)  and  Durkee  et  al  (1986)  demonstrated  this  linear  dependence  of  radiance 
on  optical  depth. 

There  are  three  primary  sources  of  satellite-detected  radiance;  aerosol  scatter  (LA). 
Rayleigh  scatter  (LR),  and  sea  surface  reflectance  (Ls).  This  relationship  can  be  expressed 
as 

L  =  LA  +  LR  +  LS  2.2 

.  Rayleigh  scatter  occurs  when  the  size  of  the  aerosol  particle  is  much  smaller  than  the 
wavelength  of  the  incident  radiation  and,  though  strongly  wavelength  dependent,  it  does 
not  vary  spatially.  Sea  surface  reflectance  is  very  small  for  red  wavelengths,  the  albedo 
is  0.5%,  and  zero  for  radiation  with  wavelengths  longer  than  0.7/tm,  except  where  the 


sun-earth-satellite  geometry  results  in  sunglint  (Ramsey  1968).  Therefore,  by  avoiding 
data  contaminated  by  sunglint.  total  radiance  can  be  approximated  by  the  sum  of  the 
radiances  due  to  Rayleigh  and  aerosol  scatter  (L  ^  LA  +  LR).  By  subtracting  the  compo- 
nent due  to  Rayleigh  scatter,  aerosol  particle  radiance  (LA)  will  remain: 

LA  =  L(cS,  n,  4>)  -  LR  ~  (  "~p  )r(&)tA,  2.3 

where  5A  is  the  aerosol  particle  optical  depth. 

With  Rayleigh  scatter  forming  a  baseline  radiance,  aerosol  scatter  is  the  dominant 
source  of  upwellcd  radiance  variations  (  AL  ~  ALA  ).  Eqn  2.3  indicates  satellite  measured 
radiance  is  directly  related  to  the  changes  in  the  aerosol  particle  characteristics  of  w0, 
p[Q)  and  t^  .  It  is  apparent  that  a  knowledge  of  p(0)  variations  is  required  to  determine 
a  more  accurate  SA  from  LA. 

B.     OPTICAL  DEPTH 

Optical  depth  is  a  measure  of  the  scattering  and  absorption  material  present  in  the 
atmosphere.  It  is  determined  by  vertically  integrating  the  extinction  coefficient  (<?„,), 
from  the  earth's  surface  to  the  top  of  the  atmosphere  (H),  expressed  by 


(5=  [    a 


xtdz.  2.4 


Extinction  is  due  to  both  absorption  {acbs)  and  scattering  {osca)  of  radiant  energy;  how- 
ever, in  the  marine  environment,  absorption  can  be  neglected,  leaving  extinction  due  to 
scattering  alone  (a„r  ^  asccr)  such  that 


ch 


(5  = 


ascatdz.  2.5 


The  amount  of  scattering  is  determined  by  the  size  of  the  particles,  their  scattering 
efficiency,  and  the  size  distribution  of  particles  in  the  column  given  by 


•oo     2 

nr  Qscai(m<  Kr)n[r)drt  2.6 

o 


where. 


nr2  =  particle  cross  sectional  area, 
Qscc,(m.  ).,  r)  =  scattering  efficiency, 

m  =  complex  index  of  refraction, 

X  -  wavelength, 

r  =  particle  radius, 
n{r)  =  number  distribution  of  particles. 

Aerosol  particles  are  most  efficient  in  scattering  radiant  energy  at  wavelengths  near 
their  radius,  becoming  less  efficient  as  the  difference  between  the  two  increases.  Any 
changes  in  the  size  distribution  will  produce  a  corresponding  change  in  SA  and  in  the 
observed  radiance  LA.  The  variations  in  radiance  between  two  wavelengths,  specifically 
red  and  near-infrared,  can  be  quantified  by  the  ratio  (Durkee  et  al.  1991) 


(LA)red         1°>0P(®)&  Aired 


5, 2  =  —r~, —  ~  ~ ,—  -, —  • 


This  ratio  (S12)  is  called  the  particle  size  index  and  is  sensitive  to  variations  in  the  slope 
of  the  particle  size  distribution,  increasing  as  the  slope  increases.  Small  particles  have  the 
greatest  influence  if  the  slope  is  steep,  while  large  particles  are  most  important  for  more 
gradual  slopes.  Generally,  (LA),ed  ,  due  to  small  particle  scattering,  is  greater  than  (LA)mr 
resulting  in  an  S]2  ratio  larger  than  1  (Mahony  1991). 

C.     WATER  VAPOR  ABSORPTION 

The  dominant  influence  on  retrieved  radiances  in  the  visible  and  near-infrared 
wavelengths  is  aerosol  scattering.  Though  water  vapor  effects  are  not  as  significant  as 
aerosol  scatter,  they  can  result  in  an  over-estimation  of  S12  which  could  be  interpreted 
as  a  change  in  aerosol  characteristics  or  on  the  other  hand,  could  mask  a  significant 
change.  As  the  desired  result  is  the  estimation  of  aerosol  particle  characteristics  from 
satellite  observations,  inclusion  of  water  vapor  effects  will  improve  the  accuracy  of  re- 
sults. 

In  Figure  1,  the  effect  of  water  vapor  can  be  seen  in  the  transmittance  (t)  for 
AVHRR  channels  1  and  2.  Note  that  channel  1  t  is  greater  than  channel  2,  indicating 
less  absorption  is  taking  place.  Therefore,  L2  would  tend  to  decrease  more  due  to  ab- 
sorption by  water  vapor  than  L,  leading  to  the  over  estimation  of  S12. 
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Figure  1.  Atmospheric  transmittance  due  to  water  vapor:  Transmittance  is 
shown  with  NOAA  &.  AV1IRR  channels  1  and  2  response  functions 
(Mahony  1991). 

D.     BR1CHTNESS  TEMPERATURE  DIFFERENCE 

An  estimate  of  the  water  vapor  content  in  an  air  column  can  be  made  using  the 
AV1IRR  10  to  13  /mi  window  region.  In  channels  4  and  5,  at  wavelengths  of  11  and  12" 
urn  respectively,  absorption  by  water  vapor  is  dominant  while  aerosol  scattering  is  rela- 
tively insignificant.  The  benefits  of  using  these  channels  is  the  colocation  and  temporal 
linkage  with  channels  1  and  2,  providing  a  water  vapor  snapshot  along  the  same  path 
as  the  aerosol  optical  depth  (SA)  measurements.  Although  not  the  ideal  sensor  windows 
for  water  vapor  content,  their  use  is  justified  within  the  accuracy  of  the  approximations 
made  (Mahony  1991). 


The  brightness  temperature  (T)  is  determined  from  the  radiance  received,  due  to 
direct  transfer,  from  the  ocean  surface  to  the  satellite  sensor.  Therefore,  the 
transmittance  (t)  is  directly  proportional  to  the  retrieved  T,  Only  absorption  reduces  the 
radiance  reaching  the  sensor;  scattering  is  relatively  insignificant. 

Dalu  et  al.  (1981)  proposed  a  method  for  determining  sea  surface  temperature  using 
brightness  temperature  difference  corrected  for  air  column  water  vapor  content.  After 
manipulation  of  Dalu's  formula,  the  vertically  integrated  water  vapor  density  (w)  is 
provided  by 

w  =  A(TA-T5)  cosd.  2.8 

where  A  is  a  constant.  Using  A  =  19600g°A'"1  nr\  a  scanning  angle  0  =  0,  and  a  variety 
of  temperature  and  humidity  profiles,  Dalu  (1986)  applied  this  relationship  to  a  radiative 
transfer  model.  The  resulting  line  graph  (Figure  2)  has  a  correlation  coefficient  of 
R  =  0.99,  with  an  error  of  ±  1.5%  m~2  in  density. 

Mahony  determined  that  the  water  vapor  influence  on  the  particle  size  index  was  a 
function  of  w  and  Sn{dry),  the  latter  being  the  magnitude  of  512  when  w  =  0.  The 
equation  for  the  curve  that  best  fit  his  data  is 

Sn{drv)  = =-,  2.9 

2  1  +  0.0332V  w 

where  Su  and  w  are  determined  from  measured  radiances. 

E.     VARIABLE  PHASE  FUNCTION 

Based  on  Frost's  work,  a  variable  single  scatter  phase  function  was  used,  having 
been  determined  to  be  more  accurate  in  its  optical  depth  estimation.  Higher  ratios,  which 
are  due  to  smaller  particles,  have  more  backscatter  than  the  lower  ratios,  due  to  larger 
particles.  This  provides  the  basic  shape  necessary  for  the  phase  function  needed  to  esti- 
mate optical  depth. 

The  two-term  Flenyey-Greenstein  scattering  phase  function  (Lenoble  1985)  was 
used. 
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figure  2.      Total  water  vapor  content  as  a  function  of  brightness  temperature  differ- 
ence (Dalu  1986). 
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where  W  is  a  weighting  factor,  £,  and  g2  are  asymmetry  curve  factors,  and  0  is  the  scatter 
angle.  All  three  of  these  factors  are  functions  of  Sn  written  as 


)V=  1.2 -(0.25  x  S12), 


2.11 


g]  =  1.04 -(0.08333  x  512), 


2.12 


g2=  1.2 -(0.58333  x  Sl2). 


2.13 


From  the  equations  2.11  through  2.13,  the  phase  function's  dependence  on  Si2  is  appar- 
ent. The  marine-type  phase  function  describes  large  particle  distribution,  low  Sl2  ratios; 
while  the  rural-type  phase  function  follows  small  particle  distribution,  high  Su  ratios. 
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As  the  ratios  change,  the  factors  change,  resulting  in  a  varying  phase  fuction  from  one 
similar  to  a  marine  phase  function  to  one  similar  to  a  rural  phase  function.  Figure  3 
shows  some  of  the  variable  phase  functions  />(0)  calculated  by  Frost  (1988)  from  the 
indicated  Sl2  values. 
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Figure  3.      Aerosol  particle  size  index  and  associated  phase  functions:      (Frost  1988). 


10 


III.     PROCEDURES 

A.     DATA  AND  PROCESSING 

The  data  examined  was  collected  by  the  NOAA  7  polar  orbiting  satellite  whose  av- 
erage altitude  was  833  km.  This  sun-synchronous  platform  crossed  the  equator  at  ap- 
proximately 0230  and  1430  local  standard  time  (north-  and  south-bound,  respectively). 
Of  the  five  spectral  bands  or  channels  on  the  AVHRR  instrument;  channels  1,  2.  4,  and 
5  were  used  in  measuring  upwelling  radiance.  Table  1  provides  the  bandwidths  for  all  the 
channels  (Jensen  1986)  plus  a  weighted  average  of  the  solar  irradiance  (Fp)  over  the 
bandwidth  used  for  channels  1  and  2  calculations. 


Table  1.     NOAA  7  AVHRR  CHANNEL  BANDWIDTHS. 

Channel 

Bandwidth 

Solar  Irradiance 

1 

0.58-0. 68  nm 

154.3536  mWcnr1 

2 

0.725-1.10  ixm 

103.0908  mlVarr2 

3 

3.55-3.93  fun 

4 

10.5-11.3  fim 

5 

1 1.5-12.5  f.im 

The  archives  of  the  National  Center  for  Atmospheric  Research  (NCAR)  in  Boulder, 
Colorado,  provided  AVHRR  Global  Area  Coverage  (GAC)  data  for  two  periods  in  April 
19S2;  the  5th  -  10th  and  the  20th  -  25th.  April  1982  data  was  analyzed  for  comparison 
with  the  work  done  earlier  by  Frost  (1988). 

All  the  data  was  processed  on  the  Interactive  Environmental  Digital  Analysis  Lab- 
oratory (IDEA  LAB)  located  at  the  Naval  Postgraduate  School,  Monterey,  CA.  The 
IDEA  LAB  produced  enhanced  images  and  graphs  from  the  processed  data. 


B.     ANALYSIS 

The  data  analysis  scheme  was  based  on  that  used  by  Pfiel  (1986)  and  Frost  (1988). 
The  water  vapor  algorithm  used  by  Dalu  (1986)  was  incorporated  to  account  for  ab- 
sorption by  water  vapor,  as  recommended  by  Mahony  (1991).  Also  accounted  for  in  the 
algorithm  were  changes  to  the  solar  irradiance  due  to  variations  in  the  sun-earth  dis- 
tance. During  preprocessing,  pixels  located  at  latitudes  greater  than  70°  or  those  over 
land  were  discarded. 

To  remove  sunglint,  an  empirical  method  using  viewing  geometry  was  applied.  This 
method,  a  weighted  sum  of  relative  zenith  and  azimuth  angles  for  the  sun  and  satellite, 
identified  approximately  the  same  contaminated  areas  as  the  Avian  sunglint  routine, 
which  is  based  on  the  Cox  and  Munk  (1954)  method. 

Ozone  absorption  effects  occur  in  the  red  and  near  infrared  wavelengths  with  the 
greatest  absorption  at  a  wavelength  of  0.69/j.m  (Fleagle  and  Businger,  1980),  which  is  in 
the  channel  1  window.  This  effect  was  inadvertently  compensated  for  twice  by  Frost.  The 
error  has  been  corrected  for  this  study. 

Each  pixel  was  then  examined  for  clouds  as  only  clear  pixels  were  to  be  processed. 

Cloud  removal  was  performed  by  the  same  method  used  by  Frost: 

1.  Test  for  bad  pixels  --  If  channel  2  or  channel  4  radiance  values  are  less  than  zero, 
after  Rayleigh  radiance  was  removed,  then  the  pixel  is  classified  as  "bad"  and  not 
processed. 

2.  Test  for  high  clouds  --  A  pixel  was  considered  to  contain  high  clouds  if  its  channel 
4  brightness  temperature  was  less  than  273°A'. 

3.  Test  for  low  clouds  -  A  pixel  was  considered  to  contain  low  clouds  if  its  channel 
2  albedo  was  greater  than  40%. 

4.  Preliminary  test  for  clear  pixels  —  A  channel  1 'channel  2  albedo  was  calculated  and 
if  that  ratio  was  greater  than  1.5,  the  pixel  was  tentatively  classified  as  clear. 

5.  Final  test  for  clear  pixel  --  The  channel  2  albedo  of  the  pixel  tentatively  identified 
as  clear  was  compared  with  the  four  pixels  above,  below,  to  the  right  and  to  the 
left.  All  the  pixels  must  fall  with  in  five  counts  of  each  other  for  the  pixel  to  be 
classified  as  clear,  otherwise  it  was  counted  as  partlv  cloudv.  A  count  is  equivalent 
to  0.1069%  albedo  (Kidwell  1986). 
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The  clear  pixels  were  then  processed  for  a  variety  of  parameters.  The  list  is  basically 

the  same  as  Frost's,  the  most  significant  difference  is  the  inclusion  of  two  parameters 

with  the  water  vapor  correction,  the  terms  include  the  word  "dry"  in  their  name. 
Channel  1 /Channel  2  ratio, 
Dry  Channel  1  /Channel  2  ratio, 
Aerosol  optical  depth, 
Dry  Aerosol  optical  depth, 
Average  phase  function  of  Channels  1  and  2, 
Channel  4/Channel  5  brightness  temperature  differences, 
Number  of  clear  pixels, 
Number  of  partially  cloudy  pixels, 
P  parameter. 

The  P  parameter  is  used  to  examine  the  phase  function  calulated  by  Frost. 

The  phase  function  p(6,S}2)  is  an  important  factor  in  the  determination  of  aerosol 
optical  depth.  Any  changes  made  in  the  phase  function  algorithm,  such  as  the  addition 
of  a  water  vapor  correction  to  S12,  will  produce  corresponding  changes  to  the  calculated 
optical  depths.  A  comparison  of  this  modified  phase  function  to  that  determined  by 
Frost  should  provide  confirmation  that  the  changes  made  produce  physically  sound  re- 
sults and  add  to  the  knowledge  of  p(6,  S,2). 

A  new  parameter  was  defined  using  the  simplified  RTF  (2.1)  that  is  calculated  from 
radiance  data  and  can  be  used  to  derive  p{6)  and  SA: 

P=   —A-  xP(0,Sl2)SA.  3.1 

ro 

Using  the  fact  that  integrating  the  phase  function  over  a  sphere  is  unity,  the  following 
form  can  be  used  to  determine  p(d,  Su): 


\p(d,Sn)ciQ  =   \Ap{e,Sn)5AdQ.  = 


APdD.  =   1,  3.2 


where  A  is  equal  to  the  inverse  of  the  optical  depth.  It  is  this  equation  from  which  the 
phase  function  will  be  calculated  for  comparison  to  Frost's  curves. 
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All  of  the  results  were  stored  and  averaged  in  1°  by  1°  boxes;  each  data  point  is  one 
box.  The  number  of  pixels  varied  from  box  to  box,  depending  on  the  number  of  partially 
cloudy  to  cloudy  pixels  found  in  each  square. 
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IV.     RESULTS 

A.  INTRODUCTION 

Ongoing  efforts  to  determine  aerosol  particle  characteristics  have  resulted  in  a 
steady  improvement  in  the  algorithm.  Water  vapor  effects  were  shown  to  be  an  impor- 
tant factor  in  aerosol  optical  depth  retrieval  but  have  not  been  previously  applied  to 
satellite  data.  In  this  study  the  latitudinal  variation  in  water  vapor  content  is  initially 
examined  to  give  an  overview  of  where  the  effects  could  be  the  most  significant.  Particle 
distribution  characteristics  are  examined  next,  both  with  and  without  the  water  vapor 
correction.  The  data  having  the  correction  will  be  refered  to  as  "dry".  Both  the  corrected 
and  uncorrected  particle  size  index  and  aerosol  optical  depth  will  be  examined  and 
compared  over  the  Pacific  Ocean  followed  by  the  latitude  profiles  of  the  averaged  values. 
Finally,  the  shape  of  the  P  parameter  curves  are  compared  to  Frost's  phase  function 
curves. 

B.  COLUMN  WATER  VAPOR 

Figures  4  and  5  show  the  average  water  vapor  content  in  an  air  column  to  provide 
a  general  overview  of  water  vapor  effects  in  the  atmosphere  for  the  periods  of  April  5  - 
10  and  20  -  25,  respectively.  The  general  structure  is  that  of  a  maximum  in  the  vicinity 
of  the  equator  with  values  then  decreasing  toward  the  poles. The  drop  is  more  rapid  in 
the  northern  hemisphere  than  the  southern  hemisphere,  the  latter  is  transitioning  from 
its  summer  season  when  temperatures  are  higher  and  more  water  vapor  is  present.  The 
greater  year-round  air  and  sea  surface  temperatures  at  the  equator  allow  more  water 
vapor  to  exist  in  a  column  of  air.  From  this  it  can  be  seen  that  the  greatest  potential  for 
error  due  to  water  vapor  is  in  the  equatorial  region.  An  examination  of  the  particle  size 
index  before  and  after  applying  the  water  vapor  correction  will  be  made  next. 
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LATITUDE  PROFILE 
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Figure  A.      Brightness  temperature  difference  for  April  5  -  10. 
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LATITUDE  PROFILE 
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Figure  5.      Brightness  temperature  difference  for  April  20  -  25. 
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C.     PARTICLE  SIZE  INDEX 

Figures  6  and  7  depict  the  uncorrected  Si2  values  across  the  Pacific  Ocean  for  April 

5  -  10  and  20  -  25,  respectively.  A  broad  region  of  high  values  extends  across  the  low 
latitudes  in  the  southern  hemisphere  in  both  figures,  indicating  the  presense  of  a  greater 
number  of  smaller  particles.  Low  values  indicating  larger  particles  can  be  seen  in  Figure 

6  across  the  northern  hemisphere.  In  Figure  7  there  are  two  noticable  maxima;  one  in 
the  vicinity  of  the  New  Hebrides  islands,  the  second  west  of  South  America.  The  north 
Atlantic  appears  to  be  a  region  of  small  particles  except  along  the  east  coast  of  the 
United  States. 

In  Figures  6  and  7,  bands  cross  each  figure  from  north-northwest  to  south- 
southeast.  These  are  a  result  of  the  moving  of  the  sensor's  aperature  to  the  left  and  right 
of  the  satellite's  path  to  view  the  side  regions.  Aerosol  scatter  is  anisotropic,  favoring 
forward  and  backward  scatter.  A  satellite  is  looking  at  backscatter  when  it  looks  directly 
at  the  earth.  But  as  the  sensor  looks  to  the  side,  it  is  no  longer  viewing  backscatter,  but 
the  smaller  scatter  associated  with  smaller  scattering  angles.  This  results  in  a  decrease 
in  the  phase  function  and  a  corresponding  increase  in  the  optical  depth.  Future  work 
will  need  to  correct  for  this  effect. 

Figures  8  and  9  show  a  comparison  of  the  particle  size  index  profiles  for  Sl2  and 
Sl2(dry)  for  early  and  late  April,  respectively.  The  SX3  profiles  in  both  figures  show  a 
broad  region  of  high  values,  approximately  2.35,  from  the  equator  to  40°S;  the  values 
then  decrease  toward  the  poles  more  rapidly  in  the  southern  hemisphere  than  in  the 
northern.  The  greatest  drop  can  be  seen  in  the  southern  hemisphere  in  Figure  8  with  a 
minimum  value  of  2.0.  The  minimum  in  Figure  9  is  2.10.  Frost's  values  were  2.0  at  the 
equator  decreasing  to  1.65  moving  southward  toward  the  pole;  a  second  smaller  decrease 
occurred  toward  the  north  pole,  dropping  to  1.9.    His  underestimation  was  due  to  a 
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Figure  6.      Particle  size  index  over  the  Pacific,  April  5-10 

computational  error  in  the  algorithm  and  over-compensating  for  ozone  absorption. 
Next,  the  corrected  Sl3(dry)  will  be  discussed. 
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Figure  7.      Particle  size  index  over  the  Pacific,  April  20  -  25 

D.     DRY  PARTICLE  SIZE  INDEX 

Figures  10  and  11  depict  the  Sa{dry)  values  across  the  Pacific  Ocean.  A  comparison 
to  the  uncorrected  S[2  values  in  Figsures  6  and  7  show  that  the  particle  size  index  values 
decreased  significantly  after  the  correction  was  applied.  Open  ocean  regions  are  now 
uniform  with  no  significant  difference  between  hemispheres. 
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Figure  8.      Comparison  of  particle  si/e  index  profiles  for  April  5  -  10:     the  solid  line 
is  Su,  and  the  dashed  line  is  Su(dry). 
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Figure  9.      Comparison  of  particle  size  index  profiles  for  April  20  -  25:     the  solid  line 
is  Sn,  and  the  dashed  line  is  Si:(Jry). 
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In  early  April,  Figure  10,  distinct  high  value  areas  remain  in  the  Gulf  of  Carpentaria 
and  along  the  coast  of  China.  Late  April  (Figure  11)  shows  two  areas  of  higher  values 
in  the  northern  Atlantic  Ocean  and  along  the  west  coast  of  North  America.  The  maxi- 
mum west  of  South  America  in  Figure  7  is  easier  to  distinguish  in  Figure  1 1  and  has 
lower  values.  The  maxima  off  the  Americas  appear  to  be  associated  with  upwelling  re- 
gions that  are  biologically  very  productive,  which  could  be  a  source  for  the  greater 
number  of  small  aerosol  particles. 

In  Figures  8  and  9,  it  is  easy  to  see  the  decrease  in  size  index  values  between  5I2  and 
Snidry).  At  the  equator  Si2(dry)  is  0.5  less  than  Sl2  in  early  April,  0.4  less  later  in  the 
month.  From  the  equator  to  the  poles,  the  difference  between  S12  and  Sl2(dry)  decreases 
to  approximately  half  the  values  at  the  equator  in  both  figures.  In  both  instances  the 
overall  shape  of  the  Sl2(dry)  curve  shows  a  minimum  near  the  equator  gradually  in- 
creasing as  you  move  toward  the  poles,  which  is  the  reverse  of  Frost's  results  where 
values  decreased  toward  the  poles.  Though  a  more  uniform  profile  was  expected,  the 
peaks  that  exist  in  the  corrected  profiles  of  both  figures  correspond  to  the  latitudes  of 
the  maxima  depicted  in  Figures  10  and  1 1. 

E.     OPTICAL  DEPTH 

The  views  of  the  Pacific  in  Figures  12  and  13,  show  the  southern  hemisphere  with 
low  optical  depths  and  the  northern  hemisphere  with  very  high  optical  depths.  The 
largest  feature  is  the  wide  area  of  high  values,  0.45  to  0.50,  that  occurs  along  the 
Mexican/Central  American  west  coast  extending  westward  into  the  Pacific.  This  feature 
is  the  result  of  the  El  Chicon  volcano  (17.30A',  95.2°  IF)  erupting  on  a  number  of  occa- 
sions during  the  period  of  28  March  to  6  April.  The  high  values  in  the  Gulf  of  Mexico 
are  also  the  result  of  El  Chicon.  The  movement  to  the  east  is  possibly  due  to 
stratospheric  winds  carrying  particles  ejected  to  those  heights  (Frost  1988).     In  early 
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Figure   10.      Dry  particle  size  index  over  the  Pacific,  April  5  -  10. 

April  (Figure  12)  the  plume  extends  only  halfway  across  the  basin,  later  in  the  month 
it  appears  to  approach  Japan  (Figure  13). 

A  second  feature  can  be  seen  over  the  East  China  Sea  in  Figure  13.  This  region  of 
high  optical  depth  corresponds  to  a  region  of  large  particles  in  Figure  7.  Gobi  desert  dust 
blowing  off  mainland  China  produced  this  feature,  as  noted  by  Frost  (1988). 
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Figure   11.      Dry  particle  size  index  over  ihe  Pacific,  April  20  -  25. 

Figures  14  and  15  show  a  comparison  between  the  optical  depth  profiles  for  SA  and 
5A{dry),  The  SA  profiles  both  show  a  maximum  value  of  0.2  at  approximately  20°N, 
which  is  due  to  the  eruption  of  the  El  Chicon  volcano.  A  second  peak  located  in  the 
high  southern  latitudes  can  also  be  seen,  The  reason  for  this  maximum  is  uncertain, 
possibly  the  result  of  sea  salt  from  the  circumpolar  winds  around  Antarctica. 
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Figure   12.      Optical  depth  over  the  Pacific  for  April  5  -  10. 

F.     DRY  OPTICAL  DEPTH 

Figures  16  and  17  depict  SA(dry)  across  the  Pacific  Ocean.  The  high  optical  depths 
across  the  northern  hemisphere  due  to  El  Chicon  and  Gobi  dust  have  greater  areal  ex- 
tent and  higher  values  than  those  in  Figures  12  and  13.  These  features  cover  regions  that 
correspond  to  the  lower  Sn{dry)  regions,  indicating  larger  particles  dominate.  Figures  14 
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Figure   13.      Optical  depth  over  the  Pacific  for  April  20  -  25. 

and  15  are  the  dry  optical  depth  profiles  dA{dry),  for  the  early  and  late  April  periods. 
They  show  that  the  water  vapor  correction  resulted  in  an  increase  in  optical  depth  val- 
ues. Increases  were  generally  on  the  order  of  0.02,  but  could  be  as  large  as  0.06,  with 
slightly  greater  changes  occurring  at  greater  distances  from  the  equator.  Significant 
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Figure   14.      Comparison  of  optical  depth  values  for  April  5  -  10:     solid  lines  arc  SA 
values,  dashed  are  SA(Jry)  values. 
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Figure   15.      Comparison  of  optical  depth  values  fur  April  20  -  25:     solid  lines  are  5/ 
values,  dashed  are  0A{ii>y)  values. 


29 


changes  occur  in  areas  where  the  number  of  aerosol  particles  are  sufficient  to  affect  the 
optical  depth. 

G.     PHASE  FUNCTION 

Figure  18  shows  the  P  parameter  as  a  function  of  scatter  angle  for  three  different 
particle  size  indices.  A  comparison  to  Frost's  plot  of  the  Henyey-Greenstein  phase 
function  (figure  3)  indicates  that  the  basic  shape  of  the  P  parameter  curves  are  similar 
to  that  shown  by  Frost.  A  smaller  particle  size  index  has  a  larger  value  in  the  forward 
direction  (between  0°  and  60°  )  than  larger  index  values.  The  larger  the  particle  size  index 
the  closer  the  curve  approaches  that  of  Rayleigh  scatter,  which  is  expected  as  smaller 
particles  dominate.  Where  the  scatter  angle  is  between  90°  to  140°,  the  P  parameter  ap- 
proaches a  minimum  value  for  each  curve.  As  the  angle  increases  to  180°  the  curves 
show  an  increase  in  P  values. 
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Figure   16.      Dry  optical  depths  over  the  Pacific  basin,  April  5-10. 
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Figure 


17. 


Dry  optical  depths  over  the  Pacific  basin,  April  20  -  25. 
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Figure   18.      P  parameter  as  a  function  of  scatter  angle:     Values  are  shown  for  three 
particle  size  indices;  1.7  (diamonds),  1.9  (triangles),  and  2.1  (x's). 
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V.     CONCLUSIONS  AND  RECOMMENDATIONS 

Frost  showed  that  his  technique  for  processing  AVHRR  data  was  effective  as  a  tool 
in  the  analysis  of  aerosol  particle  characteristics.  Though  significant  improvements  were 
made  with  his  study,  sources  of  error  have  since  been  identified  requiring  futher  im- 
provements to  his  algorithm.  The  determination  that  water  vapor  was  a  significant 
source  of  error  suggested  a  re-evaluation  of  the  April  19S2  data  in  order  to  assess  water 
vapor  impact  on  aerosol  optical  depth  retrieval. 

The  results  of  this  study  have  shown  that  inclusion  of  a  water  vapor  correction 
produces  a  change  of  approximately  15%  in  the  aerosol  particle  size  index.  Without  this 
correction,  large  regions  of  high  particle  size  index  values,  which  correspond  to  a 
predominance  of  relatively  smaller  particles,  were  actually  regions  of  lower  values 
masked  by  water  vapor  effects.  This  made  it  difficult  to  see  true  areas  of  small  particles 
as  they  were  indistinguishable  from  the  water  vapor  variations. 

Since  the  phase  function  algorithm  is  a  function  of  5I2,  the  modification  made  for 
water  vapor  will  affect  the  phase  function.  It  has  been  shown  that  the  shape  of  the 
modified  phase  function  is  similar  to  the  one  determined  by  Frost.  The  improved 
parameterization  cannot  be  directly  compared  to  Frost's  until  the  integration  (3.2)  is 
accomplished. 

The  aerosol  optical  depth  has  also  been  shown  to  be  modified  by  the  water  vapor 
correction  producing  an  approximate  change  of  5%.  The  observed  extent  and  magni- 
tude of  areas  of  high  dry  optical  depths  was  greater  than  the  optical  depths  without  the 
correction. 

Additional  improvements  to  the  algorithm  used  are  still  necessary.  Multiple  scatter 
has  not  been  fully  accounted  for.  The  assumption  used  to  allow  for  single  scatter  only 
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is  valid  for  optical  depths  less  than  about  0.2.  This  value  was  exceeded  in  the  regions  of 
high  optical  depth  associated  with  major  aerosol  particle  producing  events  such  as  the 
El  Chicon  eruptions.  Errors  in  the  values  due  to  multiple  scatter  exist  and  must  be  cor- 
rected. 

A  second  area  for  improvement  would  be  in  the  determination  of  Sl2.  Current  pro- 
cedures use  a  ratio  of  channel  1 /channel  2  radiances  to  calculate  the  particle  size  index. 
The  preferred  method  is  to  use  a  ratio  of  channel  1 /channel  2  optical  depths,  which  re- 
quires a  better  knowledge  of  the  phase  function  for  both  channels.  It  is  the  optical  depth 
which  contains  information  on  the  aerosol  particles  and  would  provide  the  best  estimate 
for  S]2.  In  this  study  the  phase  function  and  optical  depths  are  at  the  channel  1  wave- 
length. This  improvement  should  eliminate  the  banding  that  occurs  in  the  Pacific  basin 
photographs  which  possibly  hides  features  of  importance. 

Finally,  additional  studies  of  global  aerosol  particle  characteristics  are  needed.  This 
study  was  limited  to  a  single  month.  Studies  exploring  the  changes  that  occur  during  a 
year  and  those  consisting  of  a  number  of  years  by  month  or  season  would  increase  our 
understanding  of  the  role  of  aerosol  particles  in  the  global  climate. 
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